A 



ICC 2017 

WASHINGTON 
JULY 02-07, 2017 


A Novel Method for Reconstructing Broken Contour Lines 
Extracted from Scanned Topographic Maps 


Wang Feng , 1 Liu Pingzhi , 2 Yang Yun , 2 Wei Haiping 3 and AN Xiaoya 2 

1. Xi’an Research Institute of Surveying and Mapping, Xi’an, Shan’xi, China, China; bjzxgsb@163.com 

2. Xi’an Research Institute of Surveying and Mapping, Xi’an, Shan’xi, China; liupingzhi2017 @163. com, yytoall @126.com, 
axy2001 @tom.com 

3. Zhengzhou Institute of Surveying and Mapping, Xi’an, Shan’xi, China; haipingwei @163.com 


Abstract: It is known that after segmentation and morphological operations on scanned topographic maps, gaps oc- 
cur in contour lines. It is also well known that filling these gaps and reconstruction of contour lines with high accu- 
racy and completeness is not an easy problem. In this paper, a novel method is proposed dedicated in automatic or 
semiautomatic filling up caps and reconstructing broken contour lines in binary images. The key part of end points’ 
auto-matching and reconnecting is deeply discussed after introducing the procedure of reconstruction, in which 
some key algorithms and mechanisms are presented and realized, including multiple incremental backing trace to 
get weighted average direction angle of end points, the max constraint angle control mechanism based on the multi- 
ple gradient ranks, combination of weighted Euclidean distance and deviation angle to determine the optimum 
matching end point, bidirectional parabola control, etc. Lastly, experimental comparisons based on typically samples 
are complemented between proposed method and the other representative method, the results indicate that the for- 
mer holds higher accuracy and completeness, better stability and applicability. 


Keywords: scanned topographic maps, broken contour lines, reconstruction, points’ auto-matching 


1. Introduction 

Although photogrammetric interpretations of aerial photography and satellite imagery are the most common 
methods for extracting altimetry information, a great number of existing paper topographic maps, which remain the 
most frequent form of ground surface representation, provide a less expensive alternative: extracting contour lines 
from scanned topographic maps. One kind of the most interesting information shown by these paper-based topo- 
graphic maps are the contour lines, which are imaginary lines that join different points located at the same elevation. 

The aliasing and false colors are inevitable in scanned common-conditioned topographic map, a large number of 
intersects and overlaps among point, line and region features, contour lines are usually the lowest layer in topo- 
graphic map, all of those problems make the segmentation and vectorization of scanned topographic maps even 
more difficult, and also result in a great amount of caps in contour lines’ images, which could not meet the demands 
of automatic vectorization. 

At present, existing techniques dedicated to realize the reconstruction of broken contour lines still could not give 
satisfactory results yet, and the primary defects lies in the low accuracy of end points’ matching and reconnecting, 
lacking of stability and applicability for all kinds of scanned topographic maps(good, poor and worse). So the re- 
search on the reconstruction of broken contour lines is an open and interesting problem and field, it is a meaningful 
and significant work to bring forward an effective method to reconstruct broken contour lines, which could solve the 
problem of filling up the caps of broken contour lines exactly and completely. 

In this paper, we first summarize the related works in section 2. Next, in section 3, introduce the general proce- 
dure of reconstructing broken contour lines, then design a novel method to match end points accurately and automat- 
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ically, a series of key algorithms and rules for supporting this method are proposed and realized. In section 4, full 
detailed processing flow is presented. In last section, there are some comparative experiments and analysis. 


2. Related works 


In general, vectorization of contour lines from scanned map needs four major steps: ©color segmentation for 
contour lines’ image; ©binary, thinning and pruning shears; ©filling up caps in broken contour lines and recon- 
structing them; ©raster to vector conversion of all thinned and reconnected contour lines. 

Despite most of the attentions have been put on stepl and step2 for a long time, a slightly greater attention has 
been paid to step3: automatic or semi-automatic reconstruction of broken contour lines. A great amount of profes- 
sors have carried out a lot of research works about reconstruction of broken contour lines, and there are many papers 
have been published. This paper analyzed the latest related research. 

When binary image is segmentated from scanned map, Morphological filters is applied to the resulted binary im- 
age for noise reduction, and a reconnection algorithm based on local information is designed: firstly Distance is 
computed by combining Euclidean distances between extremities and differences between their directions, then ex- 
tremities in lower distance are connected (Arrighi and Soille 1999). Another approach named A* search algorithm is 
similar to Arrighi ’s method(Khotanzad and Zink 2003; Mun San et al. 2004). The images processed by those ap- 
proaches still cause some problems such as overlapping contour lines and false connections (Figure 1(a)). 

A line tracing technique for reconnecting the contour lines, in which an interruption is resolved by searching from 
an end point of a line and the corresponding end point might be found within a sector around the current direc- 
tion(Eikvril et al. 1995). 

A semi-automatic method for broken contour lines’ reconstruction based on orientation, angle and Euclidean dis- 
tance is proposed (Ismail and Sheiriff 2006). But this method is not optimal and cannot reconnect some contours 
(Figure 1(b)). 



Fig. 1 . (a) Samples of intersecting connections of contour lines(Source: Khotanzad and Zink 2003). (b) Samples of gaps unsolved 
(Source: Ismail and Sheiriff 2006) 


An approach based on Delaunay Triangulation and Voronoi Diagrams is proposed which uses the concept of me- 
dial axis of a curve X: a point in the plane belongs to the medial axis of a curve X, if at least two points located on the 
curve X are at the same distance from this point( Amenta et al.1998; Ghircoias and Brad 2001). Its disadvantage is in- 
sufficient information and unsolved gaps. 

An improved method on the basis of min pair of break-points algorithm is: colorate contour lines, build spatial 
grid-index, find connect point between two break points on the basis of different color and distance, then connect the 
connect-points(Huang et al. 2006). 

GGVF (Generalized gradient vector flow) snake model algorithm is proposed to extract contour lines and use the 
directional field information to reconnect the broken contour lines(Xin et al. 2006). An automatic approach for fill- 
ing up gaps between contour lines is: gradient orientation field generation from contour line normal; matching of 
end-points which is based on perfect matching algorithm; and filling the gap between each pair of end 
points(Pouderoux and Spinello 2007). The disadvantage of this approach is its processing time is very long for users. 

An advanced semi-automatic approach is based on local and geometric properties, such as parabolic and opposite 
directions, the y-ordinate’s differences of end points, back-tracing(Samet and Hancer 2012). But false connections 
and unsolved caps could not be avoided(Figure 2). 


ii 




U 

\\ 
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Fig .2. Samples of wrong connections and not filling caps (Source: Samet 2012) 

Knowledge based on reconnection techniques was summarized, which aim at developing various solutions for re- 
connecting broken contour lines and compare their performance based on their ability to determine the search space 
for locating the extension point and its ability to reconnect(Pradhan et al. 2013) 

A great amount of original scanned maps with high curve density are in poor or even worse quality. Reconnec- 
tion errors (intersecting and overlapping) and unsolved gaps are still not handled with effectively and completely by 
using methods discussed above. In this paper, we proposed a novel reconstruction method based on a series of key 
algorithms and rules, which aims to be a semiautomatic procedure solution for reconnecting broken contour lines 
from binary images of scanned topographic maps, avoiding wrong reconnections and increasing completeness rate 
so as to reduce the human intervention as much as possible. 


3 Proposed method 

In this section, only the task of filling up the gaps and reconstruction of broken contour lines is handled with. It is 
assumed that other necessary tasks (topographic map digitization; color image segmentation and thresholding; filter- 
ing noisy pixels; thinning and pruning the binary image) are done before the task of filling up the gaps and recon- 
nection of contour lines. In other words, an input for proposed method is post-segmentation, binary and thinned con- 
tour line image. 

The reconstruction of broken contour lines includes three stages: ©end points detection; ©end points matching; 
©reconnection of matched end points. The end points detection mean that the ending pixel has only one neighbor. 
End points matching mean that detecting and identifying the most suitable chosen end point to an unfilled end point. 
Reconnection of matched end points could be finish with Cubic Spline Interpolation or Newton Interpolation meth- 
od(Parker 1997; Giesen 2000). In this paper, Newton Interpolation method which has effective performance and 
provides users smooth curves is used. 

Based on the global topology properties(auto-enclosed, parallel and equidistant) and local properties(smoothing 
and continuous) of contour lines, in order to overcome the defects and faults of the other methods, the new proposed 
method designs and adopts a series of key algorithms, such as multiple incremental back-tracing to compute 
weighted average direction angle of end point, the max constraint angle control mechanism based on multiple gradi- 
ent ranks, combination of weighted Euclidean distance and deviation angle between end points to determine the op- 
timum matching point, bidirectional parabola control, etc. 

In order to enforce the global and local restraint for end points’ matching and reconnection, raise the accuracy 
and completion rate, and enhance the ability of dealing with all kinds of complicated samples of contour line image, 
the proposed method emphasizes its stability and adaptability throughout its procedure. 


3.1 Multiple incremental back-tracing 

Previous methods usually detect and compute end points’ directional angle by back-traceing fixed number of pix- 
els from each end point, and use the linked-line between each end point and its back-tracing terminal point to con- 
struct directional angle. Back-tracing fixed number of pixels is inadequate and may cause some problems. For ex- 
ample, the direction of end point seems to be on x-axis ory-axis. It is very difficult to match with its pair which may 
not be on the x-axis ory-ordinate due to the property that contour lines generally have opposite and parabolic direc- 
tions. In addition, the amount of back-traceing pixels is difficult to fix. If too little pixels are used for back- tracing, 
the main tendency of a contour line may be ignored (Figure 3(a)). On the contrary, if too many pixels are used, the 
information of local direction of a contour line would be difficult to be reflected (Figure 3(b)). Therefore, only one- 
time back-tracing fixed amount of pixels tends to cause wrong results and has its uncertainty. Because this, we put 
forward a algorithm of multiple incremental back-tracing and weighted mean value to get end points’ directional an- 
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gle(Figure 3 (c) (d)), aiming at applying multiple incremental back-tracing to obtain weighted mean value as end 
points’ directional angle, so that the results’ stability and accuracy are further guaranteed. 



Fig .3. (a) go back with too little pixels; (b) go back with too much pixels; (c) multiple incremental back-tracing; (d) the algorithm of 
multiple incremental back-tracing. 


In figure3(d), P x means 1st back-tracing by l x pixels, P 2 means 2nd back-tracing by l 2 pixels, P t means i-th back- 
tracing by /^pixels, P i+1 means (i+\)-th back-tracing by l i+x pixels, P n means n-th back-tracing by l n pixels. So this al- 
gorithm’s main idea is: firstly build a back-tracing pixels dataset L{l lf l 2> —,lo which stores amount of in- 

creasing back-tracing pixels, place restrictions on minimum number mm{/J and maximum number raax{/J based 
on each contour line’s actual length which end point P belongs to; next execute back-tracing to locate a terminal end 
point P/, compute linked-line’s directional angle between P and P h lastly compute weighted mean value of direc- 
tional angles’ data set f3{(3 lt (3 2 , ■■■ > Pn}> then the result would be assigned to ultimate directional angle value. 

The formulas for computing ultimate directional angle may be expressed as shown below (1)(2) (3): 

a) 

^t=i ji 

Pi = atan2((y t - y 0 )/Oi - x 0 )) (2) 

ft=W~r (3) 

Zi/=i 

Let fi represent distance’s weight of PP t , (3 t represent directional angle Zth’s back-traing, (x 0 ,y 0 ) represent start- 
ing end point P, ( x it yj represent /th’s back-tracing terminal pixel P t . 

Back-tracing pixels’ dataset L{l lf l 2 , ...,1m} should be built flexibly according to broken contour line’s ac- 

tually length. For example, set L = {3,5,7,9,11,13,15}, thus seven times’ back-tracing would be executed. Besides, 
there are some other conditions should be noticed: ram{/J = 3, max{/J = 15, with max{(p i ) < ll, ll represents 
broken contour line’s actual length. 

This algorithm takes both broken contour line’s local geometrical properties and global tendency into account. 
Therefore, its computational results possess higher accuracy, eliminating randomness and uncertainty caused by sin- 
gle back-tracing. 


3.2 Max constraint angle control mechanism 

As is shown in Figure 4(a), Let a is the directional angle of end point P(x 0 ,y 0 ), which indicates the general pro- 
gressive tendency of the broken contour line from this starting point, let /? is the inclined angle of connecting line 
between P and P'(x/,y/), which indicates the progressive tendency of PP' similarly. Thus the deviation angle 
Act, with Aa = \(3 — a\, could directly reflect the transitional smoothing degree while two matched end points con- 
necting with each other. According to the geometrical properties of contour line, matching and connecting between 
two end points ought to be gently and smoothing, the acute transitional cases, such as closed angle and acute angle, 
are disallowed (Figure 4(b) cases with A a > 60°). Therefore, it is very useful to construct max constraint angle con- 
dition to control the deviation angle, let cp represent max constraint angle, which can help to control local transition- 
al smoothing degree, and avoid abnormal transitional cases effectively, like closed angle, acute angle and so forth. 
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^ Aa < 45 
\ \ X Aa < 50° 

Aa > 60° 


Fig. 4. (a) Max constraint angle’s computation; (2) End points’ connecting and transitional cases 

The formula for computing deviation angle and max constraint angle condition may be expressed as shown below 
(4): 

Act = \p — a| < cp (4) 

The condition of max constraint angle focuses too much on single end point’s local constraint, lacking of global 
controlling function based on whole map. Maybe a given max constraint angle could solve part of broken contour 
lines’ transition when connecting them, but the caps and contour line’s shapes in the whole map are various and 
complicated, then setting single max constraint angle is deficient and not comprehensive enough. If the value is 
much bigger, this means that the constraint is too loose and liable to lead to intersecting errors. On the contrary, 
much smaller value would exclude a lot of end points which seem to be against to the condition but is fit for match- 
ing and connecting. So it is necessary to build a flexible mechanism to strength applicability, for this purpose, we 
propose max constraint angle control mechanism based on the gradient ranks. 

This mechanism’s main idea is: firstly, build a max constraint angle dataset T = {(po,(pi ,(p 2 > — ><Pn\? with 
cpi(i = 0,1,2 , ... , n) is an incremental element in T ; then use cp t as a fixed value to control transitional smoothing de- 
gree between two end points; next, move cp t to next element to control unsolved end points until the last element cp n 
has been used. 

On one side, the dataset T should be flexible and comprehensive based on considering all complicated cases of 
the whole binary image, each element’s value is ascending step by step, and restrain the smoothing degree from 
small to large. On the other side, each element cp t implies the local transitional smoothing degree under a fixed 
range. As we see in Figure 5, T could be constructed according to whole map’s morphological characteristics, for 
example, set T = (0°, 20°, 30°, 45°, 50°, 60°}, then Figure 5 shows a variety of transitional morphology under dif- 
ferent deviation angles. When cp t < 45°, the transition between two match end points is suitable and conforms to 
common cases; when 45° < (p t < 60°, the transition is unnatural and liable to cause errors, decision should be cau- 
tious when using this scope; when (p t > 60°, the transition is extremely steep, that is to say max{cpi } should be 
smaller than 60°. 

(pi>=6 0 ° 

<Pt =50° 
cpt =45° 

(Pi =30° 

(Pi = 20 ° 

0 ° 

Fig. 5. Different smoothing degree under gradient ranks’ max constraint angle 

In fact, this mechanism firstly focuses on the most gently and simple broken contour lines, and intends to accom- 
plish those end points’ matching; thereafter, it increases max constraint angle value to solve the remaining end 
points gradually. So this mechanism realizes a procedure from simple to complex, from gently to steep, reveals local 
and global geometrical properties of contour lines, could eliminate intersects and other errors. 
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3.3 Combination of weighted Euclidean distance and angle to compute matching measurement 

Previous processing algorithms have guarantee smooth transition when connecting matched end points, but there 
is still not a ultimate measurement indicator for completing right matching between either end points and reserved 
end points. In other words, it is of crucial importance to design an algorithm to provide ultimate measurement indi- 
cator. As we all know, distance and deviation angle are the two primary impact factors, whereas Samet’s meth- 
od(Samet and Hancer 2012) use the minimum y-ordinate distance as ultimate measurement indicator to decide opti- 
mal match end point. That algorithm is simple and computational-cost is very low, though, the results are unreliable. 
For example, the minimum y - ordinate distance is not equal to minimum spatial distance, conversely the minimum x- 
ordinate distance may be very large, so Samet’s method is likely to cause wrong connections (See in Figure 2 red 
drawing). Euclidean distance is the actual plane distance, and is superior to y-ordinate distance when judging near 
and far, at the same time, computing accuracy is improved dramatically. In this paper, the y-ordinate distance is re- 
placed by Euclidean distance. On this basis, we propose an algorithm which could gain a quantitative comparison 
measurement to determine optimal match end point based on weighted combination of Euclidean distance and de- 
viation angle. 

This algorithm set weightings on Euclidean distance and deviation angle partly, and combine them to compute a 
reference value, which is set as a quantitative comparison measurement, and the minimum value gained from all 
values assigns the optimal matching end point. Distance’s weighting and deviation angle’s weighting reflect its im- 
portance and influencing factor in the computing of matching measurement. The more approaching to each other, 
the more possible it is to match successfully; the smaller the deviation angle is, the more smoothing the transition of 
connection is. All in all, the minimum matching measurement has its reasonability and validity. 

The main idea of this algorithm is: computing the Euclidean distance ds t and deviation angle between the 
center end point P 0 and the chosen end points P t (i = 0,1, ... k) firstly; next computing the value D t by combining 
weighted distance and deviation angle; at last, detect the matched end points by chosing D min = min(Dj(i = 
0,1, ... /c). Finally, D t may be discretized with formula (5): 

D t = A o t * + ds t * w 2 (i = 0,1, ... k) (5) 

Let w t represents deviation angle’s weighting, represents Euclidean distance of P 0 P it , w 2 represents distance’s 
weighting. 

The measuring units of distance and angle are different, and directly combining operation is disallowed. It is es- 
sential to convert angle to distance to make them to be under the same measuring unit, so we finish this conversion. 
The main idea of the conversion is: (1) convert deviation angle to central angle of a circle under a given radius; (2) 
compute corresponding arc length; (3) replace the deviation angle with arc length to compute matching measure- 
ment. The converting formulas are below(6)(7): 

Li = \Aoi\*r = (6) 

Di = L t * w 1 + dst * w 2 (7) 

Let r represents a given circle radius, A o t represents converted central angle, L t represents arc length. 

Generally, the weightings w 1 and w 2 assigned to distance and angle should flexible and variable. Gradient ranks’ 
max constraint angle mechanism controls transitional smoothing degree effectively and excludes a lot of end points 
with acute angle, then more attention should be paid to spatial distance near or far, in case of long distance cause in- 
tersection and other errors. Therefore, distance’s weighting should be emphasized and superior to angle’s weighting, 
the latter is a beneficial supplement to the former. For example, set w 1 = 0.2, w 2 = 0.8. These two values should be 
varied according to the matching results. 


3.4 Bidirectional control based on parabola 

Although previous algorithm and mechanism focus too much on right matching and smoothing transition from 
global to local, they are insufficient to deal with end points presenting acute angle parabola shape(Figl(b) marked 
by black rectangle). After color map’s scanning and segmentting, there are a lot of end points presenting acute angle 
parabola shape in binary image. In Samet’ s method, clockwise end points’ direction codes matching algorithm was 
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designed and applied to solve this problem. But only part of these end points have been solved, a lot of end points at 
counterclockwise direction can not be solved and are liable to lead to errors. 

Bidirectional codes matching algorithm based on parabola shape is put forward in this article. Its main idea is to 
extent clockwise and counterclockwise bidirectional codes matching (see Fig6 and Table 1), in order to improve the 
accuracy and comprehensiveness in direction codes matching. This algorithm is an important supplement to solve 
complicated and acute parabola end points, raising completeness rate of whole image effectively. 

i 


2 


3 

Fig 6. Bidirectional matching of end points based on parabola shape 



Table 1 . End points’ directional code matching based on bidirectional parabolic shape 


End point direction’s code 

Selectable end point’s direction code 

1 

3 or 7 

3 

5 or 1 

5 

7 or 3 

7 

1 or 5 


3.5 The executing procedure of proposed method 

The proposed method’s general idea is that all end points are marked and fetched firstly, multiple back-tracing to 
compute one end point’s directional angle, use gradient ranks of max constrained angle to control smoothing degree 
of reconnection, and compute minimum matching measurement value to determine optimal matching end point. At 
last, Newton interpolation algorithm is used to reconnect broken contours. The detailed executing procedure flows 
are listed as follows: 

Stepl: Preprocess the color scanned map: segment and threshold color image, filter noisy pixels, thin and prun 
the binary image. 

Step2: Set threshold, initialize all parameters, including: maximum search window, extension step, weight factors 
of distance and directional angle. 

Step3: Search all end points over the whole binary image, form end point array ( P 0 , ...,P it ..., P n ], set initial max 
constraint angle (p h i = 0. 

Step4: Start the procedure of end points matching for [P 0 , . . . , P t , . . . , P n } from end point P 0 based on max con- 
straint angle cpp 

©Detect a square mask h x h and assign P 0 as the center of the square mask, go back-tracing to detect another 
end point which is on the other end of the same small contour line. If any end point comes across, do not use 
this end point, mark and exclude it. 

©Detect other end points within the h X h square mask , form a chosen end point array {P 1 ,...,Pj,...,P m }. 
©Compute P 0 ’s directional angle a, the temporal connecting line’s inclined angle (3 constructed by P 0 and its 
chosen end point P ; , Pj E [P lf ...,Pj,..., P m }, the deviation angle A o { (j = 1,2,3, ..,m) between a and P, then 
judge whether Aoj < cp t , if yes remain it, else remove it from {P 1 ,...,Pj,...,P m }, judge weather 
{P lt ...,Pj,..., P m } have been traversed, if yes form a new chosen end point array {P lt ...,Pj,..., P m r], go to 
Step @. 
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©Compute the Euclidean distance dsj between P 0 and its chosen end point Pj, Pj E { P lt ...,Pj,..., P m '}, combine 
weighted dsj and Aoy to calculate P 0 ’ s matching measurement Dj. 

©Judge whether Dj reaches to minimum. When D t reaches to minimum D min , its corresponding end point P' is 
regarded as the object point. 

©reconnect these two matched end points P 0 , P' by utilizing Newton Interpolation algorithm, and remove them 
from {P 0 , P n }, refresh the array. 

©Judge whether it remains un-filling caps in the binary image, if yes go to Step5, else go to Step 7. 

Step5: Judge whether all max restrained degree ranks have been traversed, if yes go to Step 6, else set cp k = 
(Pk+ i,go to Step 4. 

Step6: Use bidirectional parabola shape end points’ direction code matching algorithm to reconnect end points 

which possess the feature of acute angle parabola. 

Step7 : Results are summarized and programe exit. 


4 Experiments and analysis 


In order to test the effectiveness, stability and applicability, the proposed method has been applied to many dif- 
ferent samples for analyzing and comparing with other methods. Due to the paper’s length limit, three images of rel- 
atively complex topographic maps in different quality are used in this paper, and the experimental results’ evalua- 
tions are presented from two aspects: one is qualitative evaluation, the other is the analysis of measurement 
indicators, including the number of finished connections, number of false connections, completeness rate, accuracy 
rate, etc. 

We have used Samet’s method and proposed method to process three typetical images respectively. Although 
Fig7(b), Fig8(b) and Fig9(b) demonstrate that Samet’s method could implement reconnecting most of broken con- 
tour lines, obvious errors(marked with black circle), for example intersects and failing to fill up bigger or steep caps, 
are still unresolved. On the contrary, Fig7(c), Fig8(c) and Fig9(c) indicate that the proposed method has almost re- 
connected all the broken contour lines, especially realized matching almost all end points with each other correctly, 
and the little remained unconnected broken contour lines hold steep parabolas (marked with black rectangle). 




Fig. 7. Sample 1: (a) pre-processed binary image(Source: Ghircoias and Brad 201 1); (b) result of Samef s method [10]; (c) result of pro- 
posed method 



Fig. 8. Sample 2: (a) pre-processed binary map(Source: Pouderoux and Spinello 2011); (b) result of Samef s method; (c) result of pro- 
posed method 
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Fig. 9. Sample3: (a) pre-processed binary map; (b) result of Samet’s method; (c) result of proposed method 


In order to analyze and compare with other method, the following numerical indicators are used: (1) total gaps; 
(2) the number of filled gaps; (3) the number of unresolved gaps; (4) the number of wrong connections; (5) com- 
pleteness rate; and (6) accuracy rate. Here let T represents the number of total gaps, 7\ represents the number of 
filled gaps, T 2 represents the number of not filled gaps, C 1 represents the percent of reconnected number, C 2 repre- 
sents accuracy percent. In general, T is very difficult to be obtained from contour lines’ images automatically, which 
is acquired manually; T x could be recorded in processing automatically; and T 2 merely could be acquired by check- 
ing the original image and resulting image manually. The formulas of computing C x and C 2 are given below (8)(9): 


r> Ti nr | / o \ 

c i ~ |T| ( 8 ) 

r — IriMzl /q\ 
^ " I Til (9) 


Analysis and comparision results of proposed method and Samef s method are given in Table 2, and all the nu- 
merical indicators come from previous three experimental samples(Fig 8, Fig 9 and Fig 10). 


Table 2 Analysis and comparision results of proposed approach and Samet’s method 


Map name 

Map resolution 

Total gaps 

Method 

Ti 

T, -T 2 

Ci(%) 

C 2 (%) 




Samef s method 

50 

4 

94 

96 

Fig 7(a) 

581*431 

53 

proposed method 

52 

0 

98 

100 




Samef s method 

66 

3 

90 

95 

Fig 8(a) 

1278*903 

70 

proposed method 

68 

0 

97 

100 




Samef s method 

91 

6 

84 

93 

Fig 9(a) 

704*468 

101 

proposed method 

100 

0 

99 

100 


5 conclusion 

In this paper, a novel method is proposed for reconstruction of broken contour lines. Because all existed tech- 
niques and methods can still not get satisfactory results yet, the proposed method is dedicated to achieve good stabil- 
ity and effectiveness in processing different quality images. This method designs and utilizes a series of key algo- 
rithms and mechanisms based on the geometrical and topological properties of contour lines. 

The excellent contribution of the proposed method is that it can eliminate the errors of other methods, increase 
the accuracy and completeness rate in end points’ matching and reconnecting effectively. For the difficult gaps 
which hold acute angle parabola, proposed method utilizes bidirectional parabola control algorithm to solve this 
problem. At last, comparative experiments between proposed method and the other representative Samef s method 
have been accomplished, the results indicate that proposed method could get higher rate of accuracy and complete- 
ness, possess the merits of better stability, adaptability and applicability. Next work in progress of this paper is to 
decrease computational complexity and increase efficiency. 
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